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A method is proposed for studying wave and particle transport in disordered waveguide systems 
of dimension higher than unity by means of exact one-dimensionalization of the dynamic equations 
in the mode representation. As a particular case, the T — conductance of a two-dimensional 
quantum wire is calculated, which exhibits ohmic behaviour, with length-dependent conductivity, 
at any conductor length exceeding the electron quasi-classical mean free path. The unconventional 
diffusive regime of charge transport is found in the range of conductor lengths where the electrons 
are commonly considered as localized. In quantum wires with more than one conducting channel, 
each being identified with the extended waveguide mode, the inter-mode scattering is proven to 
serve as a phase-breaking mechanism that prevents interference localization without real inelasticity 
of interaction. 

PACS numbers: 71.30.+h, 72.15.Rn, 73.50.-h 



I. INTRODUCTION 

Since the original formulation of the localization problem by AndcrsonEI, the question of whether electronic states 
in disordered systems are localized at any strength of disorder or a mobility edge can be defined by a critical level of 
disorder has become a central issue. In the former case this yields an insulating-type behaviour of large samples, while 
in the latter the metallic-type transport of conducting electrons is allowed due to the existence of extended states. 

It was subsequently ascertained that the answer to this question depends substantially on the dimensionality of the 
disordered system. For the case of one-dimensional (ID) conductors in the limit of vanishing temperature one can 
prove in a mathematically rigorous way (see Ref. || and references therein) that the spectrum of the electrons, subject 
to an arbitrarily small but finite random potential in infinitely long samples, is purely discrete, i.e. all electron states 
are necessarily localized irrespective of their energy. As a result, the DC conductivity vanishes for such systems, 
whereas in finite ID samples the conductance falls exponentially as the length of the conductor grows. 

In contrast to the ID case, no mathematically rigorous theory of localization exists for two- (2DV, and tree- 
dimensional (3D) random systems. Formulation of the problem in terms of the renormalization groupu (RG) led 
to the one-parameter scaling hypothesis of localizationQ and appeared to provide a considerable progress in studying 
systems of dimensionality greater than unity. The main conviction gained from this development was that all the elec- 
tron states in both ID and 2D disordered systems are localized at arbitrary small disorder. Hence, the metal-insulator 
transition (MIT) is usually .believed to occur only in 3D systems, while (sufficiently large) ID and 2D systems are 
always Anderson insulators.u 

This opinion, established quite long ago, has xecently been challenged after unexpected experimental detection 
of MIT in dilute 2D electron and hole systems. Ej Originally obtained-on Si-MOSFETs, the results have entailed 
numerous experimental findings of MIT in other dilute 2D systemsJiJlij leading to the well-founded belief that the 
effect is of a rather general nature. To elucidate the unconventional experimental facts, different theoretical approaches 
were put forward including quite disputable ones (see the discussion in Ref. [l6]). Among the variety of options the 
greatest anticipations in explaining the experiments are mostly associated with accounting for the Coulomb interaction 
of carriers Olla Yet the transport theories for correlated particles in the presence of disorder still cannot claim for 
general acceptance because of the substantial controversy in the inteiaMeiation of the role of interaction within different 
parameter domains which correspond to diffusiveE3'E3 and localizedEirES regimes. 

Traditional understanding of electron localization in 2D disordered systems and numerous experimental facts which 
indicate the metallic nature of charge transport at low temperatures have brought about active progress in the line 
of study that is identified with quantum dephasing. The main efforts in this direction are focused on the detection of 
various phase-breaking mechanisms responsible for the delocalization of quasi-particles primarily localized by disorder. 
However, in spite of a large body of versions on offer the problem still remains open. 

Our intention in this contribution is to re-examine the conventional one-particle approach without appealing to 
scattering mechanisms other than those characteristic of quenched disorder. Evidence will be given that even within 
this elementary model the experimental findings of Refs. [HMLq, at first sight curious, are as a matter of fact quite 
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natural. Numerous attempts to interpret the results of Ref. |l0| within the framework of a single-particle approach were 
made, in particular, by improving the scaling approach.!! In this study, however, a fundamentally different strategy is 
chosen which is an alternative to the RG analysis .cil We prefer to obtain the observables directly, while conclusions 
(though indirect) about the localization of electron states are made on the basis of the results. 

It is instructive to recall that working out, even without a profound spectral analysis, practical asymptotic meth- 
ods for calculating the disorder-averaged many-particle characteristics (conductivity, density-density correlator, etc.) 
turned out to be more helpful for the establishment .of a highly advanced theory of ID random systems than the 
development of rigorous mathematical foundations.E^Ej The usefulness of such an approach can be attributed to the 
fact that in the context of the above-mentioned essentially perturbative methods one has managed to trace with the 
required accuracy the effect of mutual interference of quantum waves corresponding to multiply backscattered current 
carriers. In such a way, physical results entirely consistent with the anticipations based on mathematical predictions 
were eventually obtained. The present research was primarily induced by long-standing discontent associated with 
the lack of arguments of a comparable standard, either in favour of localization or against it, as applied to 2D systems 
of degenerate electrons subject to a static random potential. 

Commonly, the presence of inelastic scattering mechanisms is held responsible as a main cause of preventing 
quantum interference and, thus, Anderson localization.!^ Among these are the electron-phonon and electron-electron, 
interactions and other conceivable methods of energy interchange between the electron bath and the environment 
These can lead to the loss of phase (meaning energy) memory or, in other words, phase coherence of electronic states. 
Note in this connection that in the ID case the demand of energy coherence admits a large transfer of momentum 
for onefold scattering of degenerate electrons in the backward direction. This leads inevitably to considerable local 
breaking of spatial (instead of temporal) phases of the wavefunctions. It was recognized that such a large violation of 
spatial phases is quite helpful when deriving a constructive theory of ID quantum transport. On the basis of such 
arguments, the selection of efficient subsets of terms in perturbative expansion which are responsible for the effective 
interference of electronic waves scattered iteratively in a backward direction was suggested.^ This interference finally 
results in the formation of true localized states, even when one starts their perturbative construction from plain-wave- 
like trial functions which belong to a mathematical class different from that of localized functions. 

An analogous scenario of the perturbative formation of localized states in 2D systems cannot give the same result, 
since for sufficiently isotropic scattering the spatial coherence of the wavefunctions is already violated at distances 
of the order of the quasi-classical mean free path even for weak energy scattering. The coherence is maintained 
efficiently only within a small phase volume, thus .resulting in the relative smallness of the interference corrections, 
usually known as weak localization correctionstilO for non-one-dimensional systems. 

Nevertheless, the analysis of the problem of multidimensional localization with the use of a one-dimensional scenario 
turns out to be quite productive. As we shall prove below the problem of electron transport in 2D open system of 
waveguide type can be reduced without any approximations to the set of one-dimensional (though non-Hermitian) 
problems for the quantum waves propagating in individual conducting channels. The channels will be identified with 
extended waveguide modes. For the corresponding dynamic equations that are one-dimensional, an opportunity arises 
for making substantial use of spatial phases of the wavefunctions instead of their temporal parts. This turns out to be 
preferable from the technical point of view for solving stationary problems of the electron, as well as classical wave, 
transport. The reason for the usefulness of such a procedure lies in the fact that in one-dimensional problems spatial 
averaging has been shown to be highly advantageous, leading to the reduction of the perturbative expansions of the 
physical observables to a summable series. E30 

Being exactly quantum in nature, the waveguide approach used here is, to a certain degree, less obvious than 
semi-classical ones usually applied in most localization theories. Thus, in its context the clarity of the path-based 
interpretation is substantially lost. At the same time, the benefit of our method is that one-dimensional 'channel' 
equations enable one to distinguish unambiguously the coherent intra-mode scattering, which is easily interpreted 
from the standpoint of ID localization theory, and the inter-mode scattering which corresponds, although not quite 
directly, to isotropic scattering of semi-classical electrons. The quantum states in different channels are specified by 
different longitudinal momenta. This difference essentially suppresses interference of primary and scattered electronic 
waves, if they belong to different channels. As a result, the inter-mode scattering turns out to be an intrinsic origin 
of the inability for 2D electrons to be localized by weak static disorder. 

The 'one-dimensionalization' procedure suggested in this paper gives an opportunity to hiahJight the role of spatial 
coherence in the interference of electronic waves even in a single-particle approximation.!^ The results obtained 
here by conventional perturbative methods enable us to conclude about the unrealizability of the strong (Anderson) 
localization in systems whose spatial dimensionality is greater than unity, irrespective of their size. It should be 
particularly emphasized that decoherence appearing due to inter-mode scattering is unconnected to genuine inelasticity 
in the interaction of electrons with disorder. However, the difference in longitudinal energy between the conducting 
channels could be treated as a source of 'hidden inelasticity', if one is accustomed to such an interpretation. 

The paper is organized as follows. In the next section, the problem is formulated using linear response theory 
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In section |lll| we develop a method of exact one-dimensionalization which is a central point of the paper. Then, 
in section |Vj the trial Green functions supremely important for the developed technique are analyzed with the aid 
of a two-scale perturbation method. A spectral analysis of the electron system is given in section |y|. In the final 
two sections we present asymptotic expressions for the conductance and discuss the results. A pair of tedious but 
important calculations is presented in two appendices. 



II. STATEMENT OF THE PROBLEM 



A common approach used in studies of random systems of various dimensionality is to take a hyper-cube of a certain 
linear size and vary the size while searching for the conductance scaling. Such an approach seems to be natural when 
studying spectral properties of closed systems. At the same time, it is not quite appropriate for solving transport 
problems as applied to structures of waveguide type, in particular, quantum conductors of arbitrary length. 

In this paper, we examine the case of a 2D imperfect rectangular sample of length L in the x-direction and width 
D in the lateral direction y. Degenerate non-interacting spinless electrons are assumed to be confined between the 
hard-wall side boundaries y = ±D/2, whereas in the direction of current (x) we suppose the system to be open at the 
strip euds x — ±L/2. The dimensionless conductance g{L) (in units e 2 /irfi) is calculated directly from linear response 
theoryO whence at zero temperature we have 

m 4 U a 1 ' dG(r,r') dG*(r,r>) 

giL) = -TiJJ drdr — (1) 

Here the integration with respect to r = (x, y) is performed over the area occupied by the conductor 

x G (-L/2, L/2) , y e (-D/2, D/2) . (2) 

G(r, r') is the retarded Green function of the conducting electrons. Within the isotropic Fermi liquid model this 
function is governed by the equation 

[A + k% + iO - V(r)] G(r, r') = 5{r - r') . (3) 

We adopt the system of units with Ti — 2m = 1 (m is the electron effective mass), so that A is a two-dimensional 
Laplace operator, fc^ is the Fermi wavenumber of the electrons, V(r) is the 'bulk' static random potential. 

The potential V(r) in equation ([}]) will be regarded as a short-range one and not necessarily isotropic. The term 
'short-range' implies the characteristic spatial interval over which the potential is substantially varied to be small 
compared with the 'macroscopic' lengths of the problem, namely the electron mean free path and the conductor 
length. Being considered as a statistical variable, the potential V(r) will be specified by a zero mean value and the 
binary correlation function 

(V(r)V(r')) = QW{r - r') . (4) 

Here the angular brackets denote ensemble averaging and W(r) is some function normalized to unity. The explicit 
form of this function is not so important. In many cases W{r) is approximated by the delta-function, W(r) = 5(r). 
However, the method developed in this contribution permits one to consider not only isotropic and not necessarily 
local scattering events. In addition, the choice of W(r) in the form of a delta- function, apart from restricting the 
physical applicability of the model, is not quite convenient from the technical viewpoint. Thus, for example, when 



calculating the corrections to the mode spectrum, equations ( 45 ) , a formal problem can arise of the divergence of the 
evanescent mode contribution. It is certainly abaest provided the potential W(r) is not exactly point-supported — 
this problem is familiar in quantum mechanics .P 5 H 3 1 To get rid of the divergences in thc_problems of dimension more 
than one it is merely sufficient to choose the function W(r) to be less singular than i5(r).H Therefore, in order to focus 
on the main questions, in place of the correlation equality (|4|), we shall use the expression below 

(V(r)V(r')) = QW(x - x')S(y - y') . (5) 

Additionally, this form of equation (Q) allows us to consider anisotropic scattering. Similar to W{r) from (Q), the 
function W(x) in (0) is normalized to unity. 



3 



III. THE ONE-DIMENSIONALIZATION PROCEDURE 



A. The general scheme 



It seems intuitively natural for an open system with the prescribed direction of quasi-particle transport to be 
considered as to some extent a one-dimensional object. However, the realizability of 'onc-dimcnsionalization' at 
the level of dynamic equations, which would be quite important from a mathematical perspective, is not a priori 
apparent. Here the term 'one-dimensionalization' means reduction of the two-dimensional differential equation (|^) 
to a set of one-dimensional equations. Although a system of waveguide type can often be regarded as a collection of 
one-dimensional quantum channels, the latter are not independent in general. Normally, they are strongly coupled 
with each other through static or dynamic inhomogeneities present in the problem. 

Nevertheless, in what follows we intend to show that just the waveguide nature of a system under consideration 
enables the mathematical description of the transport problem for a 2D region (||) to be reduced to a set of independent 
strictly one-dimensional, although non-Hermitian, problems posed on the interval x G (— L/2,L/2), regardless of the 
strength of the disorder. To perform the reduction one should merely pass to the mode representation, i.e. Fourier 
transform in the transverse coordinate y, using some complete set {</>(y)} °f eigenfunctions of the transverse free- 
electron Hamiltonian, namely the Laplace operator in Eq. (g). The conductance (0) acquires the form 



m = 



L 2 



dxdx' 



E 

n.n' = l 



dG (x,x')dG*,(x,x') 



dx 



dx' 



(6) 



where N c = [kpD/Tr] is the number of conducting channels or, in other words, extended waveguide modes. Equation 
(|J) is then transformed into a set of coupled ordinary differential equations for the Fourier components G nn i(x,x'), 



d 2 

— + k 2 n +iO-V n (x) 



(7) 



m = l 



Here k 2 l = k\ — (nir/D) 2 is the longitudinal energy of the nth mode, U nm (x) is the inter-mode matrix element of the 
potential V(r), 



U nm {x) = / dyip n (y)V(r)tp m {y) 

J D 



(8) 



Note the difference between the summation limits in Eqs. (^|) and (R). Restriction of the summation in (|^) by 
the numb er of conducting channels implies, strictly speaking, weakness of the electron scattering, to be specified in 
Sec. [II B . Under the same assumption formula ([j]) itself is valid, where the products of the Green functions of the 
same kind (both retarded and advanced) have already been omitted. In equation (0), the summation is naturally 
performed over a complete set of waveguide modes. 

Particular attention should be drawn to the fact that the term containing the diagonal (intra-mode) matrix element 
U nn (x) = V n (x) is initially detached from the sum of Eq. ([?]), so that the matrix ||E/ nm || hereafter is held off-diagonal 
in the discrete mode variable. This little technical trick enables one to reduce the problem of finding the overall of 
the functions G nn i (x,x') to the solution of a subset of purely one-dimensional closed equations for the mode-diagonal 

functions G nn (x, x') only. To this end we first introduce the auxiliary trial Green functions G n V \x, x') (n £ N), each 
obeying the equation 



d 2 

+ kl + iO - V n {x) 



G n v \x,x') = 8{x-x') , 



(9) 



and Sommerfeld's radiation conditionsEZE3 at the strip ends x = ±L/2. These conditions seem to be natural to 
impose on an open system. For the case of the ID equation (||) they acquire the formEj 



tt T ik r 
ox 



G n (x, ; 







(10) 



x=±L/2 



It implies that the field of the nth mode radiated by the point source placed at x' G (-L/2, L/2) reaches the 
corresponding (±) end of the interval and then propagates unscattered with the conserved momentum k n beyond the 
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ends of the conductor. Possible scattering in the leads attached to the strip from the left and right should be taken 
into account separately. 

Although a solution of the stochastic problem (^|) and ( |lC)| ) cannot be obtained in quadratures, in the case of weak 
scattering specified below by inequalities ( ^p| ) the main features of the solution can be extracted by due consideration 
with any desired accuracy in the intra-mode potential V n (x) . The necessary analysis is presented in the next section. 
As for now, we merely consider all the functions G n V ^ {x, x') as a priori known ones whose properties are specified 
solely by the elastic intra-mode scattering. With these functions chosen as an initial approximation for the exact mode 
functions G ran , only the inter-mode scattering associated with the off-diagonal matrix U nm (x) will then be taken as 
a perturbation. To implement this intent, we turn from equation (0) to the consequent integral equation, 



Grin' X ) = G^(x,x')5 nn' ~r 



E 

m— 1 



dxi R nm (x,xi)G m n'{x\,x') 



Here the kernel 



R nm {x,xi) = G\f>( x ,x\) 



(ii) 



(12) 



already contains only those harmonics of the potential V(r) which are responsible for the inter-mode scattering. 
A thorough study of system ([ll]) leads to a conjecture that all non-diagonal mode elements G rnn (m ^ n) can be 
expressed only via the diagonal element G nn by means of some linear operator K, 



G mn {x,x ) — / dx\ K r 

IL 



nn 



(13) 



To specify this operator, one should separate the term with the diagonal (in the mode variable) Green function on the 
right-hand side of equation ( |ll| ) and substitute all non-diagonal Green functions in the form of equation ([l3]). Then, 
after renaming the mode variables, we arrive at the following equation for the matrix elements of the operator K: 



OO p 

l (x,x') = R mn (x,x') + ^2 / dx^mk{x,x{)K kn {xi,x') 



(14) 



Equation (H4) belongs to a class of multi-channel Lippmann-Schwinger equations that are known to be extremely 



singular, in contrast to their single-channel counterparts.Ej Nevertheless, by choosing the trial Green function G n as 
a zero approximation and perturbing it only by the inter-mode potentials, one manages to avoid the above mentioned 
singularity. Note that mode indices m and k in Eq. (|14j) take all the positive integer values except for the value n. 
This urges one to interpret the functions appearing in (|14|) as matrix elements of operators acting in two-dimensional 
mixed coordinate-mode space (x, m) which does not include the mode n (the notation M„ will be used for that space). 
The presence in Eq. (|l4|) of the right-hand index n, which does not belong to M n , can be ensured by introducing the 
projection operator P„ that will make the mode index of any operator standing next to it (both from the left or right) 
equal to n. With this convention accepted, it follows from equation (^) that the operator K implementing relation 
( |l3| ) has the form 



K 



1 — R 



RP, 



(15) 



Here R is a 2D operator acting on M„ and is specified by the matrix elements Ql^ 

As for the mathematical correctness of the operator representation (|T5|), it depends on the existence of the inverse 

operator ^1 — Rj . This point is discussed in Appendix |a], where we provide a proof that detachment of the 

intra-mode potential V n (x) in Eq. (0) prevents the possible singularity. 

Thus, equations JT3| ) and (Jl^) reduce the problem of finding the complete Green function G(r, r') within the 2D 
region (^) to calculation of its diagonal mode components only. In order to do that it is necessary to put n' — n in 
Eq. (Q) and substitute all non-diagonal components G mn from equation (|l^). As a result, the closed equation for the 
diagonal component G nn is deduced, 



5 2 
dx 2 



1 —+k 2 n + iO-V n {x)-f n 



G nn (x, x) = S(x - x') 



(16) 
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In equation ( jig ) , in addition to the local intra- mode potential V n (x) , the operator potential T n has arisen, 

% = V n U (l - R) 1 RPn = PnU (t-R) 1 P„ . (17) 



(18) 



Here U is the inter- mode operator specified on M n by the matrix elements 

(cc, to) = Ui m (x)S(x — x') . 



The potentia -3jn-,has quite an interesting interpretation. In the right-hand side of equation ( |T7j ) there is a conven- 
tional T-matrixMr 1 ! enveloped by the projective operators, one of which removes the excitation from mode n and the 
other restores it back to the same mode after the appropriate scattering events within the subspace M n . Hence, the 
potential T n , although corresponding to effectively intra- mode scattering, actually includes all inter-mode scattering 
events undergone by the excitation while passing over 'closed paths' in the mode space. The intra- and inter-mode 
scattering mechanisms turn out to be attributed to different potentials in equation (fl6|), facilitating significantly the 
subsequent interpretation of the results. In what follows the potentials V n (x) and T n will be referred to as those 
responsible for direct intra-mode and inter-mode scattering, respectively. 

Finally, in this subsection we express the conductance (^|) through the diagonal elements of the mode Green matrix. 
After rearranging the terms in Eq. (^J) and using relation ( |l3| ) we obtain 



No 



g(L) 



L 2 



n=l 



dx dx' 



dG nn (x,x') dGn n (x,x') 



m = l 



22 // dxidx2 



dx 



dx 



dx' 



G nn (xi, x )K^j n (a;, X2) 



dG* nn (x 2 ,x') 
dx' 



(19) 



Expression ( |19| ) jointly with equation (16) completes, in principle, the 'one-dimensionalization' procedure introduced 
at the beginning of this section. In this form the problem under consideration is convenient for a numerical treatment 
at any disorder strength, because the solution of the 2D problem governed by equation (||) is reduced to a finite set of 
purely ID problems @ and (|l^). At the same time, assuming weak electron scattering (in the semi-classical sense), 
we manage to proceed with our analytical consideration and obtain the results. 



B. The weak scattering approximation 

It is natural to specify the intensity of electron scattering in terms of characteristic spatial scales inherent to the 
problem. Henceforth we recognize the scattering as weak provided the following inequalities hold: 

kp 1 ,r c <l. (20) 

Here r c is the correlation radius of the potential V(r), £ = 2kp/Q is the semiclassical mean free path of electrons 
evaluated within the model of a (5-correlated 2D random potential, i.e. W(r) = S(r) in Eq. (|j). 



Estimation of the norm of the operator R specified on M n by the matrix elements (12) results in 

1,6,12 D/L 



(21) 



Under conditions ( po| ) , this enables us to find an expansion to lowest order in the impurity potential of the operator 
K, Eq. (|l5|), so that it becomes approximately equal to R almost regardless of the conductor aspect ratio. The exact 
operator T n from ([T^) can, in turn, be replaced by the approximate value 

t„ « P n UG (v ^UP n (22) 

where the operator is specified on M n by the matrix elements 

(x,k\G^ \x',m) = 6 km G^\x,x') . 



Besides the reduction of the T-operator (|17| ) to truncated form (|22j), a substitution of the approximate matrix 
elements of the operator K brings the second term of conductance (pit) to the form 



G 



j j j j i \jt { , X\) ,y\* i dG* m (x2 1 x') 
dxidx2U m n{xi)U mn {x 2 ) G y m > (x,x 2 )G nn (xi,x ) — , 



dx 



dx' 



(23) 



which is convenient for performing the ensemble averaging. It will be shown below that all Green functions in (|2£ 
and not only the trial ones, may be thought of as independent of the inter-mode potentials U mn . Yet correlation of 
those potentials with the intra- mode one, V n (x), governing the trial Green functions, can be disregarded in view of 
Eq. ([59|). As a result, after averaging conductance with the use of (||) and (|), the expression, which will be 
subject to a further analysis, takes the form 



(g(L)) 



L 2 2-~> 



dx dx' 



dG nn (x, x') dG^nix, x') 



dx 



dx' 



D ^ 



dxidx 2 W{xi - x 2 )(G i ^ ) (x,X2)^G { ^\x,xi)^(^G nn (x 1 ,x')-^G nn *(x 2 ,x') 



(24) 



IV. ANALYSIS OF THE TRIAL GREEN FUNCTIONS 

The trial Green functions Gm\x,x') enter the potential 7^, Eq. (|22|), thus determining the exact mode functions 
G nn (x, x'), and the second term of the conductance (p4|). Although these functions appear as subsidiary mathematical 
objects, they are of great concern for the problem and therefore deserve particular consideration. The study of the 
trial functions is also instructive, since it may provide useful insights into the analysis of some misinterpretation 
regarding 2D localization which existed until recently in the literature. 

Note first that the perturbative solution of equation (^|) depends substantially on whether the corresponding un- 
perturbed waveguide mode is either extended or evanescent. The Green functions of evanescent modes with n > N c 
are localized even without any perturbation. To find them in the limit of weak scattering, it is sufficient to restrict 
oneself to zero-order perturbation in the potential V n (x), 

G^)( X:X i) = -J— cxp (-\k n \\ x -x>\) , n>N c . (25) 

The problem is much more involved for the extended modes, n < N c . Inasmuch as the function G^ {x,x') is defined 
as a solution of the strictly one-dimensional problem ([}]) and (|To|) , to find it correctly in the context of localization 
theory the plain-wave-based zero approximation is not quite appropriate. This stems from the fact that in such an 
approximation it is rather difficult to account for the interference of multiply backscattered waves. Instead we apply 
the two-scale perturbation method analogous to that used in the theory of non-linear oscillationsEl This method 
showed itself well for solving-lhe problem of charge transport in extremely narrow, namely single-mode, surface- 
corrugated conducting strips.Eil Below an outline of the method is given together with some essential results. The 
details of their derivation are deferred to Appendix |B| 

The method used in Rcf. ^ is based on the representation of the ID Green function, which is the solution of the 
boundary- value problem (0) and ([l0|), via the solutions of the appropriate Cauchy problems, 

G^\x,x') = W^ 1 [tp+(x\n)ip-(x'\n)Q(x - x') + ip + (x'\n)tp^(x\n)Q(x' - x)] , (26) 

Here the functions ip±( x \n) are linearly independent solutions of the homogeneous equation (^ with the radiation 
conditions analogous to (|10| ) satisfied at the strip ends x = ±L/2, according to the 'sign' index. The Wronskian of 
those functions is W n , Q(x) is Heaviside's step function. This reformulation of a boundary- value problem in terms 
of an initial-value problem will allow one to perform averaging over the disorder later on£3 

It is advantageous to represent the functions ip±(x\n) as superpositions of modulated harmonic waves propagating 
in opposite directions along the x-axis, 

i , ±{x\n) — tt± (x\n) exp(±z/c„a;) — ij±(x\n) exp(^fik n x) . (27) 

Within the framework of the weak scattering approximation (|20|), the amplitudes 7r-|-(a;|n) and j±(x\n) vary slowly 
as compared with the 'fast' exponentials exp(±i/c„x), so that the radiation conditions for ip±(x\n) are reformulated 
as the 'initial' conditions for the smooth amplitudes as follows: 
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ir±(±L/2\n) = 1 



7 ±(±L/2|n) = 



(28) 



In view of the smoothness of ir± and 7±, differential equations for them can be derived by means of averaging the 
equations for ip±(x\n) over an arbitrary-valued spatial interval intermediate between the 'microscopic' lengths k~ x 
and r c on the one hand, and the 'macroscopic' lengths on the other hand. Among the latter lengths are the scattering 
length, to be specified in the course of the solution, and the sample length L. For weak scattering, the equations for 
7r± and 7± are reduced to the coupled first-order ones, 



TT±(x\n) ± irj n (x)w±(x\n) ± Q±( x h±( x \ n ) =° > 
K±(x\n) T ir] n {x)j±(x\n) ± ( n ±(x)ir±(x\n) = . 



(29) 



Here the variable coefficients rj n (x) and Cn±(x) are random fields associated with the intra- mode potential V n (x) in 
the following way: 



Vn(x) 



1 

1h n 



x+l 



Vn(t) , 



Cn±0*0 



1 

2k n 



x+l 



dt 
2i ( 



(30) 



x — l 



x—l 



For the intermediate property of the averaging interval 21, the fields rj n (x) and Cn±( x ) are : m fact, nothing but the 
'narrow' sets of spatial harmonics of the potential V n (x) with the momenta close to zero and ±2fc n , respectively. The 
real field rj n (x) is responsible for the 'forward' scattering, while the complex one Cn±(x) for the 'backward' scattering 
of the nth waveguide mode. 

Under the assumption of weak scattering, only binary correlators of the random potentials govern the majority 
of statistical characteristics of physical quantities. It was shown in Ref. that only the two correlation functions, 
{Vn(x)r) n (x')) and (( n ±(x)(*±(x')) , of modified random fields ( |30"| ) may be thought of as non- vanishing. Calculation 
with the use of model (||) readily gives 



(Vn(x)r]n(x')) 

(Cn±(x)C n± (x')) 



Here 



-(V) 
J f 



(n) 



2D 
3Q 



{2k n f 



(n) 



-Fi{x 



and 



2D (2k n ) 2 
3QW(2k n ) 



(31a) 
(31b) 

(32) 



are the forward and backward mode scattering lengths, respectively; W(q) is the Fourier transform of W(x) from 
Eq. (|). The function 



= ^(i-^)o(a-W) 



21 



21 



(33) 



in Eqs. (|3l]) is sharp at the scale of 'macroscopic' lengths, so it can be regarded as the ^-function in the 'distributional' 
sense, Fi(x) — > 5(x). 

Equations ( p9[ ) and correlation relations ( J3l| ) enable one to obtain the entire statistical information concerning the 
trial Green function for any extended mode n. In Appendix [b| derivation of all plain moments of that function is 
presented with the following result: 



GW{x,x') 
/i G K . 

It is noteworthy that from Eqs. 



—i 
2k n 



exp 



ijjLk n \x - x'\ - - 



3) we have the following estimate: 



-(V) 



(n). 



(34) 



-(v) 

J f,b 



(n) - NJ cos 2 i? n 



(35) 
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where d n is a 'sliding angle' of the nth mode with respect to the x-axis, | sini9„| = rni/kpO. The value of (|35|) is 
coincident, in order of magnitude, with the localization length widely believed to be characteristic for multi-mode 
quasi-one-dimensional (Q1D) waveguide systems (see, e.g., Refs. fi"4]-fl7|). However, in the present theory lengths (S3) 
are nothing but the extinction lengths of the auxiliary excitations that are not subjected to inter-mode scattering. 

Besides plain moments (p4|), important characteristics of the random function G„(x, x 1 ) are the mixed moments 



G { n\x,x') ' G ( n V> (x,x') ). At fi = v all of them are smooth (not oscillatory) functions of the argument 



2 00*, 



x- 



whose spatial decrease is determined by the one- dimensional localization length £ n = 4LP^(n)j3 The second term in 
( pi] ) contains one of the simplest correlators of this type, the 'density-current' correlator (G*VG). It was already 
studied in applications to the problem of classical wave transport .E3 We do not present here the exact expression for 
this correlator since only the fact of its exponential decrease at the localization length is of significance for our 
analysis, 

(^'(^J^^^oceip^^) ■ (36) 



V. THE MODE STATES SPECTRUM: PERTURB ATIVE TREATMENT 



For calculation of the diagonal Green function G nn (x, x') in equation ( |16|) by means of the perturbation technique, 
it would make a good sense to reconstruct the operator potential T n so that the mean value equals zero. To that 
end one has to define the result of the action of the operator (T n ) on the function G nn . For n > N c (evanescent 



modes) the function G nn can be left in its unperturbed form 
the operator nature of the potential (T n ) and reduced form 
integral has to be calculated: 



25[ ) due to the weakness of scattering. For n < N c , with 
) of the T-operator taken into account, the following 



(f n )G nn (x,x') = J dx{T n {x,xi)G nn {x u x') 



(37) 



Here the kernel T n (x, x\) is given by 



T„(x,Xi) 



E 

m = l 



U nm (x)G^\x,x 1 )U mn (xx) 



(38) 



On performing the ensemble averaging in Eq. (|38|), it is justifiable to neglect the correlation between the inter- mode 
potentials U nm {x) and the intra- mode ones, V m (x). Within the model of point-like scatterers such correlation is 
entirely absent. Yet even in the case of disorder specified by correlation function (^) it is not difficult, using the 
definition (H) and the hard-wall model of side boundaries, to ascertain the equality 

(U nm (x)V m (x')) = . (39) 

This allows us to couple the inter-mode potentials in Eq. (|38| ) only with each other, not affecting the trial functions 

G%\ x ,xi). The averaging procedure thus makes the operator (7~ n ) effectively local. 

Since equation (Jlq) is one-dimensional in the space variable, the exact Green function G n n(x, x') can be represented 



in a form similar to that used for the trial function G„ (x, x 1 ), Eq. (26). Specifically, under weak scattering conditions 
the function G nn (x,x'), being considered as a function of the first argument (x) only, is composed of two slightly 
modulated exponential summands of the appearance <fi±(x) = t±(x) exp(±ik n x). By applying the operator (T n ) to 
the functions <f>± (x) one can factor the smooth amplitudes t± (x) out of the integral thus arriving at the result 



{T n )(j>±(x) = £(&„)(/>-!- (x) , 



(40) 



where 



S (M = ^ J dx x W{x - Xi)exp [Tik n {x - x x )] (g%\x,Xi) 



(41) 



=i 
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For deriving equation (f4l|), the correlation equality was used 



(42) 



(U nm (x)U k n(xi)} = jjW(x - Xt)6 m k 

which results from definition (|J) and the correlation model (||). The sharp function W(x — x%) present in Eq. ( |41 

allows us to replace the trial functions Gm\x, x\) by the unperturbed free mode Green functions G$ (x,X\). Then, 
in the case of even W(x), the factor S(fc„) acquires the simple form 



oo 



(43) 



Here Gm (g) is the Fourier transform of the function Gm (i), and it is independent of the sign of the momentum fc„. 

From the above analysis it follows that the action of the operator (T n ) on the function G nn (x 1 x') is reduced to 
multiplying by, in general, the complex-valued quantity S(fc„). Using the explicit form of Gm\q), 



(km + lO) 2 - q 2 



(44) 



for both real and imaginary parts of the mode 'self-energy' T,(k n ) = Ak 2 — ijr^ the expressions are obtained as 
follows: 



m=l J -' 



1 



Q ^ 1 



r (v) 4D 



£ 



rfg W(q + fcn) 

2^ fc2 - Q 2 



W(k n - k m ) + W(k n + km) 



(45a) 
(45b) 



The symbol V in ( 45a ) denotes the principal value. Under the conditions of weak scattering the real part Ak\ is 
alwa ys small, |Afc^| <C k 2 , so it can be disregarded without serious consequences. At the same time, 'dissipative' term 
(451: ) plays a crucial role for the further analysis and cannot be omitted. As a result, equation ( Jlq ) takes the form 



<9 2 

^ + < + io - f„(x) - at; 



G nn (x, x') = 5(x - x') , 



(46) 



where = k^ + i/r^ an d A7^ = T n — (7^). Thus, for the analysis of equation (46) we shall regard the set of 
renormalized energies (n = 1, 2, . . .) as representing the new 'unperturbed' spectrum of the system, instead of the 
original spectrum {k 2 }. The perturbation theory can now be developed making use of the appropriate zero-mean 
potentials V n (x) and A7^ 



Note the difference in summation regions for Eqs. (45a) and (45b). The summation in (45b) turns out to be 
restricted by the number of conducting channels (extended waveguide modes), because only for n < N c the disorder- 
averaged trial Green functions in Eq. ( pl| ) are essentially complex- valued (see Eq. ([34]) at jj, = 1). The level broadening 
1 implies obligatory presence in the conductor of ot her extended modes besides the nth mode itself. In the case 



of an extremely narrow strip with A^ c = 1 the sum (45b) contains no terms, and thus the system should exhibit true 
one-dimensional properties. Specifically, the electrons in such systems can be transferred within two regimes only, 
ballistic and localized, and the conductance should go down exponentially with the length L exceeding the localization 
length 

On increasing the conductor width, as soon as the wire ceases to be single-mode (N c > 2), the situation changes 
drastically. The nth- mode spectrum acquires the level broadening (45b) and is subjected to both the potentials V n (x) 
and A7n. The physical reason for spectrum 'complexification' with the availability of more than one extended mode 
in the conductor is actually the randomization of the spatial phase of the electron wavefunction in going between the 
states with different mode energies and, consequently, mode momenta. It is just the uncertainty of those momenta 
that destroys the spatial coherence of one-dimensional quantum waves governed by equation (|46|). This destruction 
prevents interferential localization of the true mode states of conducting electrons contrary to their trial states. 
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To evaluate the 'phase-breaking' effect of the term fl45b| ) note that at any N c > 1 the estimate 1/t„ ~ Q is 
valid. In particular , in t he extreme case of a multimodc conductor (N c 3> 1) with point-like scatterers, replacing the 



summation in Eq. (45b) by an integration we obtain 

« Q/4 . (47) 

It is noteworthy that for a large number of conducting sub-bands each level width becomes independent of the mode 
number n and can thus be thought of as a universal dephasing rate inherent to the 2D conductor in g ener al. Although 



it is commonly believed that the phase breaking stems from inelastic scattering, the level broadening (45b) has nothing 
to do with such a scattering mechanism. The only important reason for the existence of an imaginary part of the 
mode energy is that the conductor has to possess more than one extended mode. 

Besides the phase breaking contribution 1 /t^ , the role of varying potentials V n (x) and AT n should also be studied. 
These zero- mean potentials can be assessed through evaluation of the corresponding Born scattering rates 1 /fn and 
I/t^K Estimation of the operator norms \\V n \\ 2 and ||A7^|| 2 yields 

Tn I J-j 



en c v ' nj. 



T, 



^min 1,— , (48a) 



%r \ min [1, — ) • (48b) 

( r ) cos 2 fin V Nd 



T, 



It hence follows that scattering caused by the non-local potential AT n is more efficient than that attributed to the 
local intra-mode potential V n (x). 

On the other hand, the inter-mode scattering caused by the operator potential l~ n in Eq. (|16|) is already taken largely 
into account through the imaginary renormalization of the mode energy in (^). From ( |48b| ) it can be seen that the 

scattering rate 1 is always small compared to the level width fl45b[ ) provided that the wire is not extremely long in 
the rr-direction, i.e. if the length L does not fall into the interval L ^> N c £. Yet even within this interval the quantity 
I/t^ cannot exceed the level broadening. If that is the case, the search for strong Anderson localization at any length 
of the multimode (N c > 2) conducting strip has no sense even without any inelastic scattering mechanisms. Probably, 
it is necessary to. invent some extra conditions to make the interferential localization possible in such conductors, e.g., 
magnetic field,E£l surface roughness,E2l etc. 



VI. CALCULATION OF THE CONDUCTANCE 



Equations (46) and (|9|) provide a way to perform a direct calculation of the disorder-averaged conductance (|24|). In 
this section, to avoid cumbersome expressions the results will be given for the case of a large number of conducting 
channels, N c ;§> 1. Nevertheless, all the estimates, as well as the final formulae, are valid for an arbitrary finite 
number of modes, N c > 1. In what follows, two summands in equation (^ ) will be considered separately inasmuch as 
the calculation techniques and the contribution of these terms in the total conductance differ significantly. The first 
summand, (g^ 1 '(L)\, will be conventionally referred to as the diagonal conductance, since it does not clearly contain 
any quantities characteristic of the inter-mode scattering. The second term in (g( 2 )(L)), will be conventionally 
referred to as the non-diagonal conductance. 

Note that the effect of phase breaking, which manifests itself strongly through complexification of the excitation 
spectrum at N c > 1, enables one to obtain the solution of equation ( p6[ ) perturbatively in the potentials V n and AT n , 
i.e. to neglect, in the leading approximation, the interference of quantum waves multiply scattered by those potentials. 
It follows from estimates ( [48| ) that such an approach is absolutely justified for conductors of length L <C N c £. Yet 
even in the case of a Q1 D w ire of length L N c £, when the potential AT n in Eq. ([46]) cannot be disregarded as 
follows from estimation ([48b|), the effect of this potential will be shown to be negligibly small. 



A. The diagonal conductance 



conductance. By applying the method of Ref. 43 



accuracy as was done diagrammatically in Ref. 25 



The presence of the local potential V n (x) in Eq. (46) does not substantially complicate the calculation of the 



this potential can be taken into account with just the same 
Accounting for this local potential leads to purely interferential 



corrections, which we are not concerned with in this paper. 
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As regards the operator potential AT ni to treat it perturbatively keeping in mind the estimate ( [48| ) it is advantageous 
to go over from the differential equation (|46| ) to the consequent integral one, 



G r , 



t (x,x') = G<®(x,x') + (G^Af n G nn )(x,x') , 



(49) 



Here Gnn(x,x') obeys equation ( [46|) with AT n = 0. In the leading approximation in the parameter L/N c £ <C 1, the 
function Gnn has the simple 'unperturbed' form 



Gnn x ) 



1 



2z ky\ 



exp 



{[ik n -l/l^]\x-x'\} 



(50) 



which nonetheless includes most of the inter-mode-scattering effects. In Eq. fl50p the mode coherence length Z„ is 



associated with the nth level broadening (45b), namely = 2k n r^f 1 . In the limit N c ^$> I its value equals 

4 v) = 4£cosi?„ . 



(51) 



Substitution of the Green function (|30]) into the first summand of Eq. ( |24| ) readily gives the diagonal part of the 
conductance 



N c J V ) 

( 5 (1) w) = E^ 



,(*>) 

<>n 

1- — exp 



L 



sinh 



,0) 



(52) 



In the limit A c 3> 1, by replacing the sum in Eq. ( p2| ) with the integral and substituting the coherence length in the 
form (p5|) , we arrive at the following asymptotic expressions for the conductance (152) : 



(g^(L))^7rN c £/L. 



(53a) 
(53b) 



In the limit of L > iV c £, it follows from ( [48q ) that the zero-order approximation in A7^ is, strictly speaking, 

(T) 

insufficient for calculating the function G n n(x, x ). Nevertheless, since the scattering rate 1/rA can only be less than 
or of the order of the level width 1 /t^ , the appropriate correction to the conductance can be reasonably estimated 
by substituting into g^(L) the function G nn (x, x') obtained from ( [49| ) to the first order in A7^. A simple but tedious 
calculation brings about the following estimation of the corresponding correction to the conductance: 



(54) 



This quantity is apparently small compared with ( 53b ). The result in (|54| ) indicates that one may account for the 
inter-mode scattering in the problem under consideration by means of the phase breaking factor which shows itself 
only in smearing of the mode energy levels. 



B. The non-diagonal term of the conductance 



When calculating the second term (</( 2 )(L)) in Eq. ( ^4|) two essentially different correlators have to be evaluated, the 
first containing the trial mode Green functions and the second composed of the exact ones. As regards the correlator 
of the trial functions Gm \ from result ( |34| ) it is clear that at L <C N c £ those functions can be replaced by the free 
ones, independently of the mode. Using the function G nn in the form (dOI) we obtain 



(9 {2 \L)) = - 



AL 2 D 



E 

71=1 



(<*>) 



&71 



• exp 



cosh ■ 



L 

7v) 



sinh 



L 

7v) 



E 



(55) 



At N c 3> 1 this expression is substantially simplified giving the asymptotics as follows: 
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L<£t ( ff P) (L) )«_^ c |, (56a) 

l<L« N c £ (g (2) {L)) w -nN c £/2L . (56b) 

In the case of extremely long conductors with L 3> N e £, the trial Green functions in (^J) cannot be replaced 
by unperturbed ones so that the effect of ID localization should be taken into account properly. The 'density- 
current' correlator standing in ( pi|) was studied in Ref. The result obtained there closely corresponds to analogous 
results for the 'density-density' and 'current-current' correlators. The main feature of the correlator is that it decays 
exponentially at the localization length £ TO = 4Lj (m). Keeping this in mind, it is not difficult to estimate the 
non-diagonal term in (^4|) as 




VII. RESULTS AND DISCUSSION 

While comparing results ([56]) and ( |57| ) with those given by equations ( ^3| ) it can be seen that the non-diagonal 
part of the conductance is not small relative to the diagonal one only if the conductor length falls within the interval 
<«L< N c £. In this case the non-diagonal term is half the size of its diagonal counterpart and has the opposite sign. 
Hence in the whole range of the conductor length (with the width being kept constant) the conductance is described 
by the following asymptotic expressions: 

(i) L < £ : (g(L)) « N c 

(ii) £<.L<^N C £ : (g(L)) » (n/2)N c £/L > 1 (58) 

(iii) N c £ <C L : (g(L)) m nN c £/L < 1 . 

Strictly speaking, result (|5^) is valid exactly when the number of channels is large, N c ^$> I. Nevertheless, even in the 
case of N c > 1 only an insignificant difference occurs, produced by the dependence of the coherence length l„ on the 
mode number n. This dependence differs somewhat in the non-semi-classical limit from that given by equation (p3l|), 
but the difference reduces to a numerical factor of the order of unity only. 

The result given by Eq. (|58| ) allows one to distinguish three regimes of charge transport in a quantum conductor, 
depending on its aspect ratio. Regime (i) corresponds to entirely ballistic transport, both from the semi-classical 
and quantum points of view. The result obtained exhibits a natural stepwise dependence of the conductance on the 
transverse size of the wire. 

In regimes (ii) and (iii) the semi-classical motion should be regarded as diffusive. This opinion is consistent with 
the conventional view of classical diffusion since the mean free path £ is small compared with the sample length L, 
although it can be in arbitrary relation to the conductor width. Furthermore, such an interpretation is supported by 
the ohmic, i.e. inversely proportional to L, dependence of the conductance in both of the indicated regimes. At the 
same time, it should be particularly emphasized that only in regime (ii), commonly called the weak localization regime, 
is the result given by the classical kinetic theory reproduced exactly. The expression for the diffusion coefficient in 
regime (iii) differs from that pertinent to regime (ii) by a factor of two. 

Regime (iii) is often called localized, because it is usually supposed that in wires of such a length (commonly 
referred to as Q1D systems) the Anderson localization should manifest itself to a considerable extent, thus leading to 
an exponential fall in the conductance. However, from the above analysis it follows that conductors with more than one 
quantum channel interconnected through scattering mechanisms, even elastic ones, should not exhibit an exponential 
dependence of kinetic coefficients on the sample length. Such a behaviour is characteristic for the single-mode wires 
only, which is entirely consistent with theoretical predictions for one-dimensional disordered systems. 

To associate findings of this paper with convictions that have prevailed hitherto, it is helpful to examine the electron 
transport in regimes (ii) and (iii) starting with the trial waveguide states governed by the homogeneous equation . 
Those states are certainly fictitious, they would exist provided the inter-mode scattering were disregarded. If so, 
the system would indeed represent a set of N c independent one-dimensional conducting channels where the true 
interferential localization should take place as a result of direct intra-mode scattering from the potentials V n . For all 



of the channels, a hierarchy of localization lengths would exist similar to that representative of equation (35). The 



length region in (ii) corresponds to the condition when the majority of the trial states are extended. In contrast, 
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in regime (iii) all of those states are localized, which is consistent with the expectation of an exponential fall of the 
conductance with the growth of the sample length. 

In reality one certainly cannot disregard the inter-mode scattering in the case of arbitrary quenched disorder. That 
scattering results in quite strong coupling of the channels or, rather, the trial mode states. The coupling has shown 
itself through the complexification of true mode spectrum in Eq. (ff6|). Such a complexification suggests that for 
any extended mode in a multi-mode strip all other extended modes can be thought of as a phase-breaking reservoir 
destroying quantum interference and hence strong (exponential) localization. Only weak localization corrections due 
to the local intra- mode potentials can be detected in the both of the diffusive regimes ( p5S| ) . A comprehensive analysis 
of the matter is beyond the scope of this paper and will be presented elsewhere. 

The existence of different diffusion regimes (ii) and (iii) can be interpreted as the dependence of the diffusion 
coefficient on the conductor aspect ratio. This dependence can hardly be extracted in the framework of the semi- 
classical approach. It results from the fact that under a gradual transition from regime (ii) to (iii), with the growth of 
the conductor length, the trial waveguide states undergo sequential localization. This should reduce the probability 
for the corresponding excitations to leave the conductor through the current terminals, whereas the probability of 
their scattering into other extended modes should increase. When all the trial states become finally localized, the 
diffusion coefficient stabilizes at the value corresponding to regime (iii) . The dimensionless conductance of such a long 
wire is less than unity as a consequence of the conventional Ohm's law. It seems that this smallness has previously 
been the reason to suppose all genuine states in Q1D conductors to be localized. 
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APPENDIX A: ON REGULARITY OF THE OPERATOR K, Eq. (15) 



To ascertain that unity is not among the characteristic numbers of the operator R, let us take advantage of the 
operator identity 

°° (-l) k 

logdet(l - R) = Trlog(H - R) = Tr V ^— ^-R* . (Al) 

fc=i 



The last equality in (Al) presumes that the operator norm is limited by ||R|| < 1. This is entirely consistent with the 
estimate (pi|). 



Consider the traces of the first two terms in sum ( Al), 



00 

TrR = / dxG^\x,x)U nn (x) , (A2) 
n=l L 

TrR 2 - jr ff dxdx'G%\ x ,x')U nm (x')G%Xx\x)U mn {x) ■ (A3) 



n,m— 1 



Since at n > N c the function G^n \x,x') in the case of weak scattering has the f orm (25), it is easy to conclude 



that the divergence of the logarithm in (Al) can arise from the first term, given in (A2), provided U nn {x) ^ 0. The 
divergence stems from locality, i.e. coincidence of the arguments, of the Green functions. It manifests itself not only 
on average but also at a given realization of the random potential entering this term. By separating the diagonal, i.e. 
intra-mode, potential V n (x) = U nn (x) and making the matrix ||?7 nm || off-diagonal we prevent the singularity of the 
operator K given by Eq. (jig). 



APPENDIX B: STATISTICAL MOMENTS OF THE TRIAL GREEN FUNCTIONS 

To perform the ensemble-averaging of the random function $^(a;, x'\n) = \Gn(x, x')] ^ , fi G K, in accordance with 
representation (p6|) and (p7|), we first decompose the function Gn (%■,%') into the sum of four terms each containing 
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narrow packets only of spatial harmonics with phases close to ±k n (x ± x r ). In doing so one should use the asymptotic 
expression for the Wronskian W n , 



W n = 2ik n [ir+(x\n)ir-(x\ri) + j + (x\n)j^(x\n)] . 

This results from the assumption of smoothness of the amplitude functions in Eq. 
and (Bl) into (26), the function G^\x, x') can be represented in the form of a scalar product 



(Bl) 



n). Then, after substituting <Y2 



Gi v) (i,x') = (e ifc " 



-ik n x 



Gi G$\ ( e 

G4 G r A 



Here Gi(x, x'\n) are the smooth amplitudes constructed from the envelopes Tt±(x\n) and j±(x\n) as follows: 

Tt-(x'\n) 



G\{x, x'\n) 
G 2 (x,x'\n) 
G 3 (x,x'\n) 
Gi{x, x'\n) 



W An{x) 



2 k n 



A n {x) 



2F An{x) 



-1 

2 k n 



A„(x) 



6- 



TT-(x\n) 

n + (x / \n) 
7r + (x|n) 

l-{x'\n) 
7r_ (x\n) 

l+(.x'\n) 
ir+(x\n) 



e _7 + y|n) r _ (a , |B) 



TT+(x\n) 

e+r + (x\n) 



7_(x'|n) 
7r_ (x\n) 



TT + (x'\n) 

B_ ; 1 ; I -\X\n) 



n + (x\n) 



The notation used in Eqs. (B3) is 



A n (x) = [l + r + (x\n)T-{x\n)Y 



r±(x|n) 



6+r+(x|n) 



l±{x\n) 
7r± (x\n) 



7r_ {x'\n) 
7r_ (x\n) 



G± = Q[±(x - x')] 



(B2) 

(B3a) 
(B3b) 
(B3c) 
(B3d) 



Before averaging the functions ( p33|) over the random potential, let us note some useful features of the dynamic 
system (p^). Since Tr±(x\n) and j±(x\n) are the causal functionals of the fields r] n (x), Cn±{x) and (^(x), they 
are determined by the values of those fields on the intervals (x,L/2\ and [— L/2, x) for the functionals labelled by 
the indexes (+) and (— ), correspondingly. The Green function elements ( |B3[ ), which will be subjected to ensemble 
averaging, are constructed in such a fashion that supports of the random functions entering the functionals of 'plus' 
and 'minus' type do not meet. Due to the random fields being effectively ((-correlated, see Eqs. ( |3l|) and (|33|), averaging 
of the functionals with different si gn indexes can be performed independently. 

It also follows from equations (29) and conditions (Eq) that all terms of functional series for n±(x\n) contain an 
equal number of fields C«± and Cn±-> whereas "f±(x\n) has an extra functional factor Cn±- Since for weak scattering 
all fields r) n (x), Cn±(x) and £*_j_(a;) are approximately Gaussian random processes, only the first summands in square 

brackets of Eqs. (B3a) and ( B3b|) r emain non-zero after averaging, while the quantities vanish. By the same 

arguments, the factor A n (x) in (B3) can be replaced by unity. 

In view of the statistical independence of the functions r] n (x) and Cn±(x), it is convenient to average over the real 
field rj n (x) already at the initial stage. To that end, it is advantageous to perform the following phase transformation 
of the amplitudes tt± and r y±: 



±L/2 

Tj n (xi) dxi 

±L/2 

rj n (xi) dxi 



TT±(x\n) — w±(x\n) exp ±i 
l±{x\n) = 7±(x|n) exp =p i 

The new amplitudes tt± and 7± obey the equations 

V ± {x\n) ± Q±(x)j±(x\n) = 
7± (x\n)± Cn± (x)tt± (x\n) = 

where the random field £ n ± ( x ) is related to ( n ± (x) by the equality 



(B4) 



(B5) 
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Cn± (x) = Cn± (x) exp 



•±L/2 

± 2i I r] n (xi) dxi 



(B6) 



This latter condition does not modify correlation properties of the backscattering fields, Eqs. (J3l|) . Then performing 
a Fourier transformation of the function $ M (a;, x'\n) over x' , we arrive at the expression conveniently decomposed into 
the sum of 'plus' and 'minus' functionals, 



<$>^(x,q\n) 



2k r 



,,iqx 



<f>l + \x,q\n) + $( ] {x,q\n) 



Here the functions <f>^ are given by 

,±L/2 

$&>(x,q\n) = ± / dxi 

J X 



TT±(xi\n) 



tt± {x\n) 



exp 



-iq(x — xi) + ifxk n \x — x\ | ± ifi / r] n (x2) dxi 

J xi 



Averaging functions (JBq) over the random field rj n {x) with the use of (31a) readily yields 



,±L/2 

&lt'(x, q\n) — ± / dx 

J X 



7r±(a;i|n) 



TT±(x\n) 



exp < —iq(x — x%) + 



\x-x\\ 



(B7) 



(B8) 



(B9) 



To then perform averaging over the fields C, n ±{x) it is convenient to use the dynamic equations for the functions 
^^\x, q\n) and T±(x\n) — j±(x\n)/ir±(x\n). They read 



d$f\x,q\n) _ 1 
dx 



2L { p{n) 



$V>(x, q\n) - /iCiWri^ln)^)^, g|n) 



(BIO) 



± 



dT±{x\n) 
dx 



<n±{x) + (*±(x)T 2 ±(x\n) 



These equations stem from definitions (B8) and system (B5) along with the obvious 'initial' conditions 

¥±\±L/2, q\n) = 0, T±(±L/2\n) = . 

Averaging of (BIO) with the use of Furutsu-Novikov formulaS gives the equation 



dffi\x,q\n)) 



dx 



= 1 - 



1 



2 \Lf\n) 



- ifik n =F iq 



(Bll) 



(B12) 



(B13) 



from which the result (p3) arises immediately. 
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